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We derive lower bounds for the attainable fidelity of standard entanglement purification protocols 
when local operations and measurements are subjected to errors. We introduce an error parameter 
which measures the distance between the ideal completely positive map describing a purification 
step and the one in the presence of errors. We derive non-linear maps for a lower bound of the 
fidelity at each purification step in terms of this parameter. 



I. INTRODUCTION 

Entanglement purification |l]-|3| is one of the most important tools in the theory of Quantum Information and, in 
particular, in Quantum Communication. It allows, in principle, to create maximally entangled states of particles at 
different locations, even if the channel that connects those locations is noisy These entangled particles can be 
then used for faithful teleportation Q or secure quantum cryptography 

The basic idea in entanglement purification is to "distill" a few N' pairs of particles (qubits, for example) in highly 
entangled states out of N > N' pairs in a mixed state with lower fidelity of the entanglement (or, in short, fidelity) 
using local operations and measurements. This fidelity is defined as the maximum overlap of the density operator 
of a pair of qubits with a maximal entangled state. If the initial pairs are in a non-separable state [||J^], then one 
can obtain asymptotically (in the limit N — > oo) maximally entangled states [Io| provided all local operations and 
measurements are perfect [Q,[nj. In practice, there will be errors both in the local operations and measurements. 
The purpose of this paper is to analyze this problem for the purification protocols introduced in Refs. HQ. We are 
interested in analyzing the conditions under which one can purify in the presence of errors, as well as in the limitations 
of the purification protocols. In particular, we find a non-linear map which relates a lower bound for the fidelity at 
two consecutive steps of the purification protocol, which allows us to derive lower bounds for the reachable fidelity. 
In order to analyze this problem, we introduce a parameter <5 which characterizes the errors. It measures the distance 
between the ideal operations and measurements and the ones in the presence of errors. 

Quantum Communication in the presence of errors has been previously considered by Knill and Laflamme jl2| in a 
general context, and by Van Enk et al. |l3f| for a particular experimental setup |14| , The work of Knill and Laflamme 
introduced ideas of fault-tolerant quantum computation Jl5| to show that there exists an accuracy threshold for 
storage of quantum information, which also applies to the case of Quantum Communication. As shown by Bennett 
et al Q one can rephrase this result in terms of entanglement purification with one-way classical communication. 
In Ref. |16j| , entanglement purification together with a generic error model is used to estimate the possibilities of 
quantum communication over long distances using quantum repeaters. The employed entanglement purification 
protocols explicitely utilize two-way classical communication, which makes them much more efficient for quantum 
communication. In the present paper we use purification protocols which utilze two-way classical communication, 
and therefore our lower bounds are much higher than those derived from the theory of Knill and Laflamme fll2| . On 
the other hand, we are interested in a rigorous lower bound for the achievable fidelity for arbitrary errors, and not 
in an estimation |ll|. The results and methods developed here can be generalized to derive lower bounds for other 
interesting problems in which local operations and measurements are imperfect, such as quantum teleportation or 
quantum cryptography. 

This paper is organized as follows: Section II contains a summary of the main results of this paper, and is directed 
to the reader who is interested neither in the technical details of the definitions of our error parameter, nor in the 
derivations of the non-linear maps for the lower bound of the fidelity. In Section III we introduce the error parameter 
5 and derive some properties related to the fact that it is a distance between completely positive linear maps. Finally, 
in Section IV we derive the non-linear map for the fidelity of entanglement in terms of this distance and sketch its 
dynamics. 
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II. SUMMARY OF THE MAIN RESULTS AND DISCUSSION 



In the standard scenario of entanglement purification M, two partners at different locations share N pairs of qubits, 
each pair being in a state described by a density operator p. A purification procedure produces N' < N pairs in a 
state p' "closer" to a maximally entangled state ^> me by only using local operations, local measurements, and classical 
communication between the partners. More specifically, if we define the fidelity of the entanglement 

F(p) = max(V'me|p|V'mo), (1) 

where the maximization is taken with respect to maximally entangled states ij) m ei then F(p') > F(p). In the following 
we will call F(p) simply fidelity. 

It has been shown p0[ | that if p is non-separable (it cannot be written as a convex combination of factorized density 
operators |^||) then there are purification procedures which obtain F(p') = 1 in the asymptotic limit N — > oo. In 
particular, if F(p) > 1/2 one can reach this goal by using the purification procedure devised by Bennett et al Q] 
and improved by Deutsch et al 0. It consists of a concatenation of purification steps involving two pairs of qubits, 
which give rise to a single pair with higher fidelity. In all these procedures, one assumes that the local operations and 
measurements are error free. In a real situation, however, there will be errors due to the coupling to the environment, 
imprecise apparatus, etc. Although small, they will limit the maximum attainable fidelity and will dictate whether 
purification is possible or not. 

In this section we first briefly review the purification protocol introduced in Refs. 0,0], and define the notation 
that we will use later on. Then we consider the same procedure in the presence of general errors, and characterize 
these errors in terms of a single parameter, d, which basically expresses the departure of the purification step in the 
presence of errors from the ideal one. Next, we express the lowest possible fidelity (worst case) in each purification 
step as a function of the lowest possible fidelity in the previouos step, which leads to a non-linear map. We analyze 
this map and discuss the conditions required for purification with imperfect means. The properties of our definitions 
and the technical details are presented in the following sections. 



A. Error free purification protocols 



In this subsection we review the two purification procedures introduced in Refs. |lL|?j. Subsequently we will refer to 
them as Scheme I and II, respectively. We characterize them in two different ways: firstly, in terms of a completely 
positive linear map between the initial density operator and the one after the measurement; secondly, in terms of 
a non-linear map relating the diagonal matrix elements of the density operator at each step in the Bell basis with 
the ones in the previous step. In the next subsection we will generalize the first characterization to the case of 
imperfect operations in order to introduce the parameter describing the errors, and then we will generalize the second 
characterization to find a lower bound for the fidelity. 

The purification protocols I and II both consist of a sequence of steps in which local operations are applied to two 
pairs of qubits, followed by a measurement of one of the pairs which is then discarded. Depending on the outcome of 
the measurement, the other pair is discarded or not. In the latter case the fidelity F\ of the remaining pair is larger 
than that of the original ones. This step is applied to the N pairs obtaining N-y < N/2 pairs of higher fidelity F-y. 
Then it is applied to the resulting iVi pairs obtaining N 2 pairs of fidelity F% > F\. Continuing in this vein, one can 
reach asymptotically F n — > 1 when n — > 00. 

Let us consider a single purification step. It starts out with two pairs 1 and 2 in the state p\ 2 — P <S> p, applies the 
local operations described by the superoperator IA (unitary in the case of Scheme II) and then measures each of the 
qubits of the second pair in the basis {|0), |1)}. We denote by x the outcome of the measurement: x = if the qubits 
are found in the state |0) 2 = |00) 2 ; x = 1 if they are in |1) 2 = |H)2i x = 2 if they are in |2) 2 = |01) 2 ; and x = 3 if 
they are in |3) = |10) 2 (the subscript 2 denotes the second pair). We denote by V x (x = 0, . . . , 3) the map defined as 
follows 



V x (p 12 )= 2 (x\U(p 12 )\x) 2 . (2) 

This map is linear and completely positive. The probability of obtaining the outcome x is p x {p\ 2 ) = tr [V x (pi 2 )]- If 

the outcome is x = 2, 3, then the first pair is discarded and otherwise it is kept. In the latter case, the state of the 
first pair will be 

, Vo(P12)+Vi(P12) ... 

Pi = : : : 7-. (3) 

PoKPW +Pl(Pl2) 
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Thus, each step of the purification protocol is completely characterized by the maps Po,i- (Note that V x stand for 
different maps depending on whether we are discussing Scheme I or Scheme II.) 

On the other hand, if one is only interested in the fidelity at each step, one can use a simpler characterization 
of each purification step in terms of four real numbers. In the purification protocols I and II, the local operations 
characterized by U consist of a bilateral XOR gate and specific single qubit rotations. In that case, the diagonal 
elements of the density operator p' in the Bell basis only depend on the diagonal elements of the density operator 
p, and therefore each purification step can be characterized by a non-linear map between these four diagonal matrix 
elements. We denote by A l n = {(f) l \p n \(j) 1 ) where p n is the density operator after the n-th purification step and \<j> 1 ) 
are the elements of the Bell basis (i — 0, 1, 2, 3), 

|^ 4 ) = -J=(|00)±|11)) 

|0 2 > 3 ) = -±=(|OI)±|IO». 

In particular, = F ni the entanglement fidelity at each step. For Scheme II there is according to Ref. a simple 
non-linear map that relates A n +\ to A n , namely 



(P\V (Pn ® Pn) + Vl{ Pn <g> Pn )W) 

(4) 



tr [Va(p n ® Pn) + Vl(pn ® Pn)} 
f(A n ) 



g{A n ) 



where 



f°(A n ) = (AW + (Al l )i, (5a) 

f 1 (A n ) = 2AlAl (5b) 

f(A n ) = (Air + (Alf, (5c) 

f(A n ) = 2A° n Al (5d) 

giA^^iAl + Alf + iAl + Alf. (5e) 

The map (0) has a fixed point at A = (1, 0, 0, 0), which is reached if the initial state has Aq = F > 1/2 [T^] . This fact 
expresses that in the absence of errors, one can use this purification protocol to purify states with F > 1/2 and reach 
a fidelity as close to one as we please. 

Scheme II |]J is governed by a similar map. The main difference is that at the end of each step the resulting state 
is brought into Werner form, that is the three diagonal elements A 1 , A 2 , A 3 are made equal to (I — A°)/3. Therefore 
one can concentrate on the first diagonal element, the fidelity A , only. The fidelity after the nth purification step is 
then given by 



f°(A° n 



l-A" 



A - K 3 ; (P.) 
n+1 ~ 9{aV-^)' (6) 

As (§, this map has an attractive fixed point at A = 1, and all Aq > 1/2 are attracted to it. 



B. Characterization of errors 



In practice, while performing the purification protocols errors will occur, in the local operation as well as in the 
measurements. The imperfections in the local operations can be accounted for by substituting the action of the 
superoperator U in Eq. (^) by the action of some other completely positive, trace preserving linear map. The errors in 
the measurements will be related to the following fact: in practice, the outcomes x = 0, 1 will be ultimately attributed 
to the presence/absence of clicks in some kind of detectors. Due to imperfections, the projection operators (or, more 
generally, POVMs) corresponding to these clicks are not exactly the same as the ideal ones (||). Consequently, the 
probabilities of the outcomes x = 0, 1 as well as the state remaining after the measurement will differ from the ideal 
ones. In general, we can describe both these erroneous operations and measurements in terms of a single completely 
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positive linear map V x which does not necessarily preserve the trace (we will use tildes in the case in which there 
are errors to distinguish them for the error free case). That is, if the two pairs are initially in the state pn = p® p, 
a purification step yields the outcome x with a probability p x {pi2) = tr[P x (pi2)]. The state of the pair after the 
measurement is 

7Wi 2 )+7Mpi 2 ) f7 . 
Pi = -=—, — s —/ — r- (') 

Pq(P12) +Pl(Pl2) 

Thus, as before, the maps Vo,i completely characterize each purification step. 
We characterize the errors by a single parameter as follows: 

5 := mzxd(V x ,V x ), (8) 

x=0,l 

where d{V,V) denotes a distance between V and V . The explicit form of this distance is given in Eq. ( [l3|) below. 
We emphasize that for a given set-up, one can (in principle) perform local measurements to completely characterize 
"P x , and therefore obtain the value of 8 experimentally The error parameter 8 has a clear physical meaning 

since it measures the distance between the ideal process and the erroneous one. We would like to remark here that 
due to the fact that there are measurements and postselection involved in the process, we have to work with maps 
V x that do not preserve the trace. In Section III we discuss why it is adventageous to use those maps instead of trace 
preserving maps. 

Some remarks concerning the adopted description of errors are in order: We envision V as the reduced dynamics 
of the two entangled pairs coupled to some environment. As shown in p(| reduced dynamics need in general not be 
completely positive (not even positive) on the whole system space. In taking the imperfect system dynamics to be 
completely positive we do (as discussed in pfl|) essentially assume that there is no initial entanglement between the 
system and any environment to which it might be coupled during gate operations. There may be, however, initial 
entanglement of the system with another environment that is not affected by the gate operations. As in the error-free 
purification schemes [yj^] we also assume the two pairs that participate in a purification step to be disentangled from 
each other. 



C. Purification with imperfect means 



Once we have defined a parameter that characterizes the errors at each purification step, we can analyze the whole 
purification procedure |y,[7J in the non ideal case. In order to do that, we define A l n — (</> l |/5„|</> 1 ) where p n is the 
density operator after the n-th purification step. We are particularly interested in the fidelity at each step A° = F n . 
In Section IV we show that for suitable initial conditions Aq and error parameter 8 

A° n >a n , Ai<b n , (n = l,2,...) (9) 

where 

a 2 n + bj-2S 

a ' n+1 (on + K? + (1 - a n - b n f + 26 liUaJ 

K+1 al + (l- a n f - 26 (1 ° b) 



and aa — Aq, bo — Aq. For Scheme I only the fidelity A° and therefore the bound (10a) with b n replaced by (1 — a„)/3 
is relevant. 

Equations ([Toj ) define a non-linear map that can be iterated to yield a lower bound for the attainable fidelity 
-Poo > Ooo which depends on the value of 5. In the following we will analyze the map (|Tc|). 

Let us first concentrate on the fixed points (a/, 6/) of this map, and consider in particular Scheme II. In Fig. I 
(solid line) we have plotted a/ as a function of the error parameter 6. For small values of 5 < 0.01 there are three 
fixed points. The ones with largest and the smallest value of a/ are attractive, whereas the intermediate one is a 
saddle point attractive in one direction and repulsive in the others. For larger values of (5, only the smallest one 
survives. This means that for the appropriate initial values of a and b$ il 5 < 0.01 one increases the fidelity using 
the purification protocol II to a value larger than the one given by the right wing of the appropriate curve of Fig. 1. 
For example, for 8 ~ 0.005 one can obtain a fidelity F > 0.95. 
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Now, let us analyze for which initial conditions (ao,&o) the map converges to the fixed point with the largest a/, 
i.e., for which purification is possible. In Fig. 2 we have plotted in the (a, b) parameter space the curve (separatrix) 
between the stable regions for several values of 5 (5k — 0.002fc, k = 0, 1, . . . , 5). For any initial value (do, bo) lying to 
the right of each curve, the map will converge to the corresponding fixed point (asterisks in the plot). For S = 0.006 
(k = 3 in the plot), for example, one can purify from values of ao > 0.69 up to values of F > a,f < 0.94; for 5 = 0.002, 
one can reach F < 0.98 starting from ao ^ 0.61. The results show that the error threshold for purification is much 
less restrictive than the one for quantum computation |L2| . 



III. DISTANCE BETWEEN TWO POSITIVE MAPS 



We denote by H a finite dimensional complex Hilbert space and by L(H) the complex Banach space of linear 
operators A : H — > H with the trace norm \\A\\ = trQA^ A\ 1/2 ) = tr(\A\) (as usual, \A\ = {A^Al 1 ^ 2 ). We denote by 
C(H) C L(H) the convex set of positive linear operators p acting on H with ||p|| < 1, and by P(H,H') the set of 
completely positive linear maps V : C(H) — > C(H') fulfilling 

ITOII<IHI (11) 



For positive operators, the trace norm simply coincides with the trace, and therefore Eq. (11) is equivalent to 

tv[V(p)]<tt(p)<l. (12) 

Given two completely positive maps V ,V £ P(H, H'), wc define their distance 

d(V,V)= max\\V(p)-V(p)\\. (13) 
peC(H) 

It is straightforward to show that d is indeed a distance by using the fact that the trace norm is a norm. 

With this definition, we can characterize the errors by using the parameter 5 as defined in (||) . The motivation for 
this definition with respect to other possible definitions is that it easily gives lower bounds even for physical processes 
where there are measurements and post selection (as it is in the case of entanglement purification, cf. next section), 
i.e. when the map describing the physical process is not trace preserving. On the other hand (although we will not use 
this property here) it allows to easily bound the distance between processes which are composed of several individual 
processes in terms of the distances between the individual processes themselves (see next subsection). 

One can define other distances between trace preserving maps: for example, one can consider the map V that 
transforms pn — ► p' Xl where p[ is given in (Q) in terms of the linear maps Vo,i- This new map, although trace 
preserving, is nonlinear. If one defines distances between V' and the corresponding (trace-preserving) ideal map V' , 
problems related to the non-linearity arise: for example, it can happen that while the distance S between the linear 
maps T 5 , V is very small, the similarly defined distance between the non-linear maps 7-", V is of the order of 1, which 
makes the definition useless to derive bounds. The reason is that low probability processes get "magnified" by the 
normalization and then dominate the maximization used to define the distance. 

One can still define other error parameters to find sharper bounds to the fidelity in entanglement purification. 
However, by increasing the number of parameters one does not gain too much and the bounds become more complicated 
to analyze. On the other hand, d(V (8 l,V (8 1) ^ d(V, V) [ ^9[ , which would allow us to use d in processes for which 
the system in which we perform operations and measurements is entangled with another system, without having to 
include the other system in the error analysis. This may be useful, for example, in quantum computation where 
operations are performed on single qubits that are entangled with many other qubits. In that case, one can define 
other distances, as it is done in Ref. E3]. In any case, in quantum communication if we can bound the fidelity when 
the system is not entangled, we can automatically derive a bound for the entanglement fidelity [fl2l|J . 



A. Properties of d 

In this subsection we derive some properties of the distance d introduced above. Given V ,V £ P(H,H') we have: 
(1) We can restrict the maximization in (|l3| ) to one dimensional projectors, i.e. 

d(v,p)= max \\vmm-vmm\- (14) 

V>G-H\|||^)|| = l 
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Proof: We just have to prove that the distance as given in ( |14| ) is always larger or equal than the one given 
in (|l3|), since the converse is clearly true. For any p £ C(H) we write p = ^2Pi\fa)(fa\ with J^-Pj < 1 and 
ipi normalized states of H. Using the linearity of V and V and that || -Pi-AjH — max ill^M|j we nn d that 
1 1 V (p) — V (p) 1 1 < maxi 1 1 V ( | fa ) (fa | ) — V ( 1 fa ) (fa \ ) \ \ ■ Taking the maximum with respect to p in this inequality completes 
the proof. □ 



(2) For all p £ C(H) and <j> £ H (normalized state) we have 

(faV( P )\fa-d(V,V) < (faV( P )\fa< (faV( P )\fa+d(V,V) 
tr [V(p)\ - d(V, V) < tr \v(p)} < tr [V(p)\ + d(V, V) 



Proof: For (15a) we use 



\(faV(p)-V( P )\fa\ < \\V(p)-V(p)\\ < d(V,V), 



whereas for (g_5bj) we use 



tr 



V(p)-V(p) 



< tr 



V(p)-V(p) =d(V,P). 



(15a) 
(15b) 

(16) 

(17) 
□ 



Next, we give a property that allows one to bound the distance when one applies sequential maps. This may be 
useful when one has a concatenation of processes. 

(3) Given V £ P(H', H") and Q £ P(H, H'), we define V o Q £ P(H, H") according to (V o Q)(p) = V[Q(p)]. Then, 
we have 



d(V o Q, V o Q) < d(V, V) + d(Q, Q). 
Proof: Using the properties of a distance, we have 

d(V oQ,VoQ)<d(VoQ,VoQ)+d(VoQ,Vo Q), 

On the one hand, we have 

d(VoQ,VoQ)= max \\V[Q(p)] - V[Q(fi)]\\ 
pec(H) 

< max \\V(p')}-V(p')\\=d(V,V) ) 
p'eC(H') 

where we have used (|ll|) for Q. On the other hand, 



d(VoQ,VoQ)= max \\V[Q(p)] - V[Q(fi)]\\ 

P ec(H) 

= max \\V[Q(p)-Q(p)}\\. 



(18) 
(19) 
(20) 



(21) 



Now, since Q(p) — Q(p) is self-adjoint, we can substitute in this last equation its spectral decompostion 

Q( P ) - Q( P ) =X>>M Mfi(p) - Q{ P )\4>) (22) 



obtaining 



dCP oQ,PoQ)= max V 
pec(H)^ 

< max y 



\Q(p)-Q(pM WW) 
\Q(p) - Q(p)\4>) 



(23) 



max \\Q(p)-Q( P )\\=d(Q,Q), 

p£C(H) 



(24) 



G 



which completes the proof. □ 

(4) Finally, we show that the distance d stems from a norm, which may be useful to derive some other properties. First, 
let us enlarge the set C(H) so that it becomes a Banach space. The simplest way is to define S(H) — \iiifi{C(H)}, that 
is, the set of operators that can be written as a (finite) linear combination of positive operators with real coefficients. 
The real Banach space S(H) C L(H) is simply the space of self-adjoint operators acting on H . In the same way, 
we can enlarge the set P(H,H'). First, given a map V S P(H,H') we define V : S(H) — * S(H) by using the 
linearity of V [that is, if S(H) 3 A = J2i^iPi with Pi € C(H), we define V(A) = E< V(Pi)]- Then > we define 
Q(H,H') — lmn{P(H, H')} which is a real vector space. Using the operator norm 

||P|| p = max \\V(A)\\, (25) 

P A£S(H)\\A\\<1 

it becomes a real Banach space. With this definitions we have 

d(V,V) = \\V-V\\ ap . (26) 



Proof: We show that the distance given in (|13|) is smaller or equal than the one defined in ( |13| ) , since the converse is 
obviously true since C(H) c S(H). For any A € S(H) with < 1 we can write A = J2i^i\<ft}(4>\, where J2i M = 1 - 
Now, arguing as in the proof of the property (1), we obtain that ||P(A) — P(A)|| < max^ | |"P(|0) (<f>\) — V(\(f))(^>\)\\. 
Taking the maximum over all possible A £ S(H) we complete the proof. □ 
The distance d is not unrelated to other quantities used in the literature to characterize erroneous operations. 
Typically, given one of the other quantities, one can bound d (and vice versa within the respective domains of 
applicability) . Specifically this is true for the minimum fidelity, the error amplitude ]l2] | , and the generic error model 
p6| . The diamond norm introduced in [ p)[ is a generalization of the distance used here and particularly useful to 
discuss operations on systems that are strongly entangled with other systems. 



IV. NON-LINEAR MAP FOR ENTANGLEMENT PURIFICATION 



In this section we derive the non-linear map (10) for the bounds of the diagonal matrix elements in the Bell basis 
of the density operator after each step of the purification process. Let us denote by A % n = (<fi l \p n \4> t ),i = 0..3, where 
p n is the density operator of a pair of qubits at the n-th step. Analogous to (||), we have 



Using (FTa) we have that 



((f> l \V (Pn ® pn) + Vl{p n ® p n )\P 



1 n+l 



tr 



Vo(Pn ® Pn) + V\{p n ® Pn) 



(27) 



f(A n )-2d 



< K+x < 



25 



(28) 



where /* and g are defined in 
in detail. 



g(A n ) + 26 g{A n ) - 25 

In the following subsections we will discuss the two purification schemes separately 



A. Scheme I 

As stated above for the scheme I we can use Eq. (^ instead of /° and forget about the other three diagonal elements. 
This gives 

A o > ( i ") 2 + ( 1 #) 2 -^ (29) 

{Al + '-^Y + (1 - Al + 1^)2 + 25 

Now we observe that the rhs of (^9|) is monotonically increasing with Jl" for all A^ > 1/8. Therefore replacing ^4" by 
| < a-n < A^n in ( p9| ) yields a lower bound for A n+1 . Since the interval [1/8, 1] is mapped into itself by the lhs of p9|) 
we arrive at the dynamical system defined by ag — A^ and 
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0>n+l 



+ (1^)2 _ 25 



K + + (i 



25' 



(30) 



For every n the value of a n is a lower bound of the fidelity after n purification steps. 

In the case 5 = the original map of Bennett et al is recovered. The three fixed points of that map at ai(5) 
0.25, ai(<5) ss 0.5, and a u (5) « 1 survive even for nonzero 5 and are given by the roots of the cubic polynomial 

x 3 - -x 2 
4 

They are plotted as a function of 5 in Fig. 1 (broken line). For 5 > 0.008 only the lower fixpoint survives. 

The upper and lower fixpoints are attractive, while the intermediate is repulsive. Consequently even an imperfectly 
implemented Scheme I allows to purify ensembles with initial fidelity Fi„ > di(5) up to a fidelity F out > a u (5), 
provided that 5 < 0.008. 



"7 


9 / 




"1 








x — 






8 " 


"1 S _ 




8 





B. Scheme II 



Scheme II converges faster than Scheme I and can tolerate somewhat larger errors, but the analysis becomes 
significantly more complicated, since all four diagonal elements of the density matrix come into play. Using we 
have 



A > 



{AD 2 + (Al) 2 26 



(A° n 



< 



2A 2 n A\ + 26 
{Al+AlY + iAl+AlY-25- 



(31a) 
(31b) 



To proceed the same way as in the previous subsection we need again a monotonicity property of the right hand 
sides (rhs) of the Eqs. (31) so that we can replace the values A l n (which are typically not known, since their exact 
value depends on the unkown errors in V) by lower or upper bounds, resp. 



Using Y^i A % n = 1 we can express the rhs of (31a) in terms of A^A^ only. It is straight forward to check that the 
resulting expression is monotonically increasing in A^ and monotonically decreasing in A\ for all , A^) fulfilling 



Al > 



■36 



1 - 25 



and Ai < 0.5. 



(32) 



Thus provided that A^ > a n , A\ < b n , and (a„, b n ) fulfill the condition ( |32| ) then a n+ i as given in Eq. (10a) is a lower 
bound for A^ +1 . 



It remains to justify Eq. ( 101 ) . Starting from ( [31b| ) we can this time express the rhs only in terms of a n = A 2 n + A^ 
and p n = A^ — A^ using the normalization condition: 



A 1 < 



a 2 n + (l-a n ) 2 -25' 



Now it is easy to check that the rhs of this inequality is monotonically increasing in a n (for fixed (3 n ) and takes (for 
fixed a n ) its maximum at j3 n = 0, where we use the fact that a n < 1 — A° and A° > 0.5. Since a n = A 2 t + A^ < 
1 — A® n < 1 — a n we arrive at Eq. (10b) by replacing f3 n — > and a n —> 1 — a n . 

The discrete dynamical system defined by the map ( |Io|) has for < 5 < 0.01 three fixpoints with a-coordinate 
around a; « 0.5, ps 0.6, a u w 1. Figure 1 (solid line) shows them as a function of 5. For 5 > 0.01 only the lower 
fixpoint survives. The exact a values are given by the real roots of a polynomial of seventh degree or equivalently by 
the intersections of the curves b n+ i(a) and 



bfix(a) = -a + \ a - (1 + 



2a- 1 



(33) 



the latter of which is defined by a n +i(a n ,bfi X (a n )) = a n . The corresponding 6-coordinates are b n +i(a x ), where 
x = l,i, u. 



8 



As in the previous case the upper and lower fixpoints are attractive, while the intermediate one is now a saddle 
point, attractive in one direction and repulsive in the others. Now essentially the same argument as in the previous 
subsection applies: points between intermediate and upper fixed points are purified to a final fidelity F out > a u . There 
are, however, two complications: First, the eventual fate of a point (a, b) depends on both a and b. Second, we need 



to make sure that the conditions (32) arc fulfilled in every step of the iteration, otherwise it is no longer valid to 
interpret (a n ,b n ) as bounds of the actual values (A°,A*). For both of these complications we have been unable to 
find complete analytical answers. Therefore we first give the numerical results before mentioning partial analytical 
solutions. 

Numerical calculations show that the physically meaningful set {(a, b) : < 1,0 < b < 1 — a} is divided in two 
parts by a curve passing through the intermediate fixed point, the separatrix (see Fig. 2). Points to the right of that 
curve converge to the upper fixed point, points to the left towards the lower one. Moreover the points to the right 
do all satisfy the conditions ( |32] ) and so do the orbits of all these points. For all ensembles described by density 
matrices with diagonal elements Aq,^4J in that region, a n ,b n as defined in ( fTo"| ) provide lower and upper bounds for 
the respective fidelities after n purification steps. For initial values to the left of the separatrix our approach allows 
no statement. The case S = in Fig. 2 indicates, how many "good" points our worst-case-consideration misses: as 
shown in Jl7t the exact border of the set of purifyable points in the (a, 6)-plane is given by the straight line a — 0.5. 

For a subset of the points to the right of the separatrix it is easy to prove convergence: All the points (a, b) fulfilling 
a > cii,b < bi, and a + b < 1 converge to the upper fixed point P u (except for Pi, of course). 

Proof: The proof proceeds in four steps. The main tool is the monotonic dependence of a ra+ i,6„ + i on a and b. 
(It is easily checked by calculation, that the coordinates of the intermediate fixed point satisfy the conditions (|3^) for 
all S so that monotonicity holds.) 

(i) Consider (a, b) in the set enclosed by the two curves b n+ i(a) and bf ix (a) (Ep. |33|, cf. Fig. 3). For these points, 
we have for all n 

a n +i > a n and b n+ i < b n . 

Since a n and b n are bounded by the coordinates of the upper and intermediate fixpoints, they form monotonical, 
bounded sequences and converge therefore. Since a n increases and b n decreases, they converge towards (a u ,b u ). 

(ii) Similarly it is seen that all points (a > a u , b < b u ) do converge to the fixed point "from above". 

(iii) Now, consider a point A = (a, b < b u ) below the curve &„+i(a). 

Let us call a point (a, 6) better than (a',b'), if a > a' and b < b'. Monotonicity implies that if (a, b) better than 
(a',b') then this will also be true for the images of these points after one iteration of the dynamical system. 

Now compare A with X' = (a' = a, 6') between the curves but with the same a as A, and with X" — (a" > a u , b" = 
b). Clearly, A is better than A' but worse than A". Since both A' and A" converge towards the upper fixpoint, so 
does A. 

(iv) A similar argument applies, if we compare a point Y = (a, b > bfi X (a)) with Y' — [a! < a.b' = b) between the 
curves and Y" = (a" = a, b" < b) below the curves: the primed points converge to the upper fixpoint, and thus (a, b) 
- being better than Y' and worse than Y" - does so, too. This completes the proof. □ 



V. SUMMARY 

The entanglement purification protocols 0,0 in the presence of errors in gate operations and measurements have 
been investigated. The errors are quantified by a single parameter derived from the trace norm. We have shown that 
these protocols allow to increase the fidelity of the entanglement even if implemented with imperfect quantum gates 
and measurements, as long as the errors are below a threshold of the order 1%. We derived a non-linear map to 
calculated a lower bound for the fidelity after n purification steps. A polynomial is given, the root of which gives a 
lower bound for the asymptotically attainable fidelity. 

This work was supported in part by the Osterreichischer Fonds zur Forderung der wissenschaftlichen Forschung 
and by the European TMR network ERB-FMRX-CT96-0087. G.G. thanks Wolfgang Diir for useful discussions. 
Part of this work was completed during the 1998 Elsag-Bailey - I.S.I. Foundation research meeting on quantum 
computation. 
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FIG. 1. The fixed points of the non-linear map: the in- 
tersections of a horizontal line at 5 with the plotted curve 
give the a-coordinates of the fixed points for Scheme I 
(broken) and Scheme II (solid). 



8=0.0 




8=0.002 
8=0.004 
8=0.006 
8=0.008 



8=0.01 



FIG. 2. The solid lines show the border between the 
two stable sets (the separatrix) for six values of 8. The 
asterisks show the corresponding (S increasing from right 
to left) upper fixed points. 
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